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Abstract. The reionization of the universe by stars and 
quasars is expected to be a highly inhomogeneous process. 
Moreover, the fluctuations of the matter density held also 
lead to inhomogeneities of the free electron distribution. 
These patterns gave rise to secondary CMB anisotropics 
through Thomson scattering of photons by free electrons. 
In this article we present an analytic model, based on our 
previous work which tackled the reionization history of 
the universe, which allows us to describe the generation of 
these secondary CMB anisotropies. We take into account 
the "patchy pattern" of reionization (HII bubbles), the 
cross-correlations of these ionized regions, the small-scale 
fluctuations of the matter density held and the contribu- 
tion from collapsed objects. 

For an open universe, we hnd that the angular cor- 
relation function C(0) displays a very slow decline from 
C(0) ~6x 10~ 13 up to the scale 9 ~ 10~ 3 rad where 
it shows a sharp drop. On the other hand, the power- 
spectrum 1(1 + l)Ci/(2n) (and the "local average" Si) ex- 
hibits a plateau of height ~ 10~ 13 in the range 10 3 < 
I < 10 6 . We hnd that for large wavenumbers I > 10 4 the 
signal is dominated by the contribution from collapsed 
halos while for I < 10 4 it is governed by the large-scale 
correlations of HII bubbles. This implies that one cannot 
discriminate reionization by stars from a quasar-driven 
scenario since the size of ionized regions never dominates 
the behaviour of the anisotropies. Moreover, the secondary 
CMB anisotropies arise from a broad range of redshifts 
(7.5 < z < 10 for the IGM and < z < 7 for galac- 
tic halos). Thus, we hnd that the generation of these 
anisotropies involves several intricate processes and they 
are close to the resolution limit of current numerical sim- 
ulations. The signal expected in our model might bias 
the cosmological parameter estimation from CMB exper- 
iments such as Planck or MAP, and could be detected by 
future mm- wavelength interferometers (e.g., ALMA). 
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1. Introduction 

Observations of the spectra of distant quasars show that 
the universe is highly ionized by z — 5, while recombi- 
nation took place at z ~ 1100. In current cosmological 
scenarios, the reionization (and reheating) of the universe 
occurs at z r ; ~ 10 (typically 6 S5 z r i ^> 15) when struc- 
ture formation is sufficiently advanced to build a large 
number of radiation sources (galaxies or quasars) which 
photoionize the IGM (e.g., Valageas & Silk 1999a). How- 
ever, the whole reionization history is a gradual and in- 
homogeneous process: each emitting object builds an HII 
region in its surroundings and reionization occurs when 
these bubbles overlap. This last stage is very rapid (e.g., 
Gnedin 2000) but at earlier redshifts there is a very inho- 
mogeneous phase which evolves rather slowly, as the size 
of the ionized regions grows and the number of radiation 
sources increase. Then, this process can leave an imprint 
on the CMB through Thomson scattering of photons from 
free electrons. First, the mixing of photons corning from 
different lines of sight leads to a damping of small-scale 
primary fluctuations. Second, the Doppler effect (photons 
get some of the peculiar momentum of free electrons) gen- 
erates secondary anisotropies since the distribution of free 
electrons is highly inhomogeneous. Thus, observations of 
CMB anisotropies could provide some information on the 
properties of the reionization process and on the features 
of the IGM at high redshifts. 

As pointed out by Sunyaev (1978) and Kaiser (1984) 
the oscillations of the velocity field (as opposite sides of 
overdensities have almost opposite velocities) lead to a 
strong suppression of these secondary anisotropies. How- 
ever, the modulation produced by the spatial variation of 
the number density of free electrons removes this cance- 
lation on small scales and can generate signiheant CMB 
anisotropies. These fluctuations of the density of free elec- 
trons can be produced by several processes. First, as ex- 
plained above, spatial variations of the ionized fraction of 
hydrogen due to patchy reionization provide a source of 
inhomogeneities (even if the IGM is uniform) . This is rel- 
evant before reionization. Second, the fluctuations of the 
matter density held itself lead to inhomogeneities of the 



P. Valageas, A. Balbi & J. Silk: Secondary CMB anisotropies from the kinetic SZ effect 



density of free electrons. This occurs both before and after 
reionization. When the density fluctuations are in the lin- 
ear regime this corresponds to the Ostriker-Vishniac effect 
(Ostriker & Vishniac 1986) while the non-linear regime is 
usually called the kinetic Sunyaev-Zel'dovich effect (e.g., 
Sunyaev & Zel'dovich 1980). 

In this article, we study both processes (patchy reion- 
ization and matter density fluctuations) in a unified fash- 
ion. To this order, we use an analytic model described in a 
previous paper (Valageas & Silk 1999a) which we built to 
investigate the reionization and reheating history of the 
universe. It includes a model for galaxy formation (de- 
scribed in details in Valageas & Schaeffer 1999) and for 
the quasar multiplicity function, which have been com- 
pared with observations at low redshifts (z < 4.5). More- 
over, it also provides a description of the correlations of the 
matter density field which is consistent with these mass 
functions. The underlying model of the non-linear density 
field is based on the stable-clustering ansatz as detailed in 
Balian & Schaeffer (1989) (see also Valageas & Schaeffer 
1997). This allows us to take into account density fluctu- 
ations within the IGM, the reionization process through 
the creation of HII regions and the correlations of these 
ionized bubbles. In addition to the IGM, we also consider 
the contribution from galactic halos. Here we restrict our- 
selves to the temperature anisotropies and we do not con- 
sider polarization. Thus, the main goals of this article are 
to: 

- present an analytic model which can describe in a 
more detailed fashion than previous works the generation 
of these secondary CMB anisotropies. Moreover, our ap- 
proach is self-consistent and it agrees with observations at 
low redshifts. 

- obtain the redshift distribution of the contributions 
to these anisotropies. 

- take advantage of the fact that we use an analytic 
model to explicitly separate the contributions from dif- 
ferent physical processes (small-scale density fluctuations, 
patchy reionization, correlations of ionized bubbles). This 
allows us to see which information can be recovered from 
observations of the CMB. 

- investigate whether one can discriminate reioniza- 
tion by stars versus quasars from the observed CMB 
anisotropies. 

- check whether the main features of this process 
strongly depend on the cosmological scenario. 

We present in Sect.|| the formalism we use to com- 
pute these secondary CMB anisotropies, as well as the 
approximations we introduce. Then, in Sect.|3j we describe 
our numerical results for the case of an open universe. 
This corresponds to the model we used in a previous work 
(Valageas & Silk 1999a) to investigate the reionization his- 
tory of the universe. In the last section we also consider 
briefly a critical density universe for comparison. 



2. Secondary anisotropies 

2.1. Contribution from inhomogeneities of the free 
electron distribution 

As described for instance in Gruzinov & Hu (1998) and 
Knox et al. (1998), Thomson scattering of CMB photons 
off moving free electrons in the IGM generates secondary 
anisotropies. Thus, for small optical depths the tempera- 
ture perturbation At (7) = ST/T on the direction 7 on 
the sky is: 



A T (7) = - 



dl.— o~Tn e fe 
c 



arc 



dx 7-V 



1 



-n e ,i 



(1) 



where ut is the Thomson cross-section, n e> { the number 
density of free electrons, v the peculiar velocity and I the 
coordinate along the line of sight, all in physical units. 
We also defined x as the dimensionless radial comoving 
coordinate: 



dx 



Az 



y/n A + (1 - n m - n A )(i + z) 2 + n m {i + z y 



(2) 



Since in our scenario reionization occurs at z — 6.8, the 
optical depth t from x — up to \ in (fil) is very small 
(t < 0.024, see Fig. 15 in Valageas & Silk 1999a) and it 
plays no role so that we used the approximation e~ r ~ 
1. Then, the two-point correlation function C(9) of these 
temperature distortions is simply: 



C{6) 



0~T C 



dxi 



dX2 



1 + Zi J l + z 2 



„ / 7l-Vl 72 - v 2 > 

x ( n e ,f,i n e ro) 



(3) 



where the directions ii j/ 1 and 7 2 make the angle 9. One 
can distinguish two effects which contribute to the cor- 
relation function in (|3|). Firstly, a uniform reionization 
in a homogeneous IGM (i.e. n e { 1 shows no fluctuations) 
provides a non-zero contribution through the fluctuations 
((7i- v i)(72- v 2)) of the velocity field. However, because of 
the oscillations of the velocity correlation (related to the 
infall from opposite sides into potential wells) the inte- 
gration along the line of sight leads to a strong suppres- 
sion of these Doppler effects (Sunyaev 1978; Kaiser 1984). 
Secondly, the inhomogeneities of the free electron number 
density add a second contribution which avoids this cance- 
lation of the velocity term. Most of the power which builds 
the fluctuations of the velocity comes from scales of the 
order of i?_i where the local slope of the power-spectrum 
is n = — 1. This gives R-i ~ 10 comoving Mpc. On the 
other hand, the characteristic scale of the inhomogeneities 
of the electron distribution is of order i?„ e ~ 100 comov- 
ing kpc (it is set by the size of the ionized bubbles and the 
scale of non-linear structures). Thus, we make the approx- 
imation that the velocity fluctuations are not correlated 
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to the density field and we write: 



C{6) 



d X idx 2 (l + ^i) 2 (l + 2 2 ) 2 ( 



7l-Vl 72 -V2, 



x [((1 + 5i)» e ,i(l + 5 2 )x e>2 ) - ((1 + <5i)a; e ,i)((l + S 2 )x e , 2 }} 



(4) 



where x e = n e j/n e is the ionization fraction, S = (rib ~ 
rib)/nb is the density contrast and we defined: 



To = 






(5) 



Here n e ^ is the present number density of free electrons: 



n b po ( A y 

n e ,o = tz 1 - — 



p 



(6) 



where Y = 0.26 is the helium mass fraction. The dif- 
ference which enters the term in the second line of ex- 
pression (EJ) corresponds to the fact that here we only 
consider the contribution to secondary CMB anisotropies 
due to the modulation of the velocity field by the fluc- 
tuations of the free electron number density. Thus, for a 
uniform free electron density field the quantity C(9) writ- 
ten in (||) vanishes. As we discussed above, the factoriza- 
tion of the velocity and density averages we used in (EJ) 
is valid at small scales r <C i?_i where we actually have: 
((7i-vi)(7 2 -v 2 )) - ((7i-vi)(7 2 -vi)) ^ ((7i-vi) 2 ). To get 
the last approximation we used the fact that the lines of 
sight ■j 1 and 7 2 must be nearly parallel since all scales of 
interest are much smaller than the Hubble scale. In other 
words, at small distances we can consider the velocity field 
to be constant. However, in practice we shall also con- 
sider larger scales where the cross-correlation of ionized 
bubbles generates some power. These scales are close to 
R-i hence we keep the expression (H) rather than making 
the approximation ((7i-vi)(7 2 .v 2 )) ~ ( v ' 2 )/3- However, 
on these scales we can expect to underestimate somewhat 
the angular correlation C(ff) since the velocity field should 
be correlated with the density field. We note £, v (Re + 1) 
the velocity term and £ e (Re + 1) the electron density cor- 
relation which enter the expression (Eh and we write (0) 
as: 



C(0) 



/TT POO 

d x -f(l+zf J d/&(R e +lKe(R0+l).(7) 



Here I is the physical length along the line of sight while 
Rg is the physical distance at redshift z between both lines 
of sight (Re is orthogonal to 1) : 



Rg{z) 



c V(z) 
#o T+z 



where the dimensionless comoving angular distance T> is 
defined by: 



V(z) = 



1 



v 7 ! - a; 



ft. 



sinh I y/l — fi m — 17a 



X 



(9) 



In this article we only consider the secondary anisotropies 
due to the inhomogeneitics of the free electron number dis- 
tribution, described by (R). As discussed above, the con- 
tribution which arises from the fluctuations of the velocity 
field only is smaller because of the cancelations along the 
line of sight. 

2.2. Velocity fluctuations 

First, the velocity term in (|7]) is obtained as follows. Since 
most of the scales which contribute to the velocity fluctu- 
ations are in the linear regime until z = 0, we can use the 
linear relation between the velocity and density fluctua- 
tions (Peebles 1980): 

v(x) = aH(z)f[n m (z)} [ d 3 fc e lkx % 6(k) (10) 



k 2 



where /(f2 m ) ~ f2^ 6 if J7a = 0, a is the scale factor (a = 1 
at z = 0), k and x are comoving quantities. Since we 
consider lines of sight which are almost parallel (small 
angles 8) we write the velocity correlation as (e.g., Groth 
et al. 1989): 

A POO 

&(x) - — (aH(z)f[n m (z)]f / dfc P{k- z) 
c Jo 

ji(kx) 



jo{kx) 



kx 



(11) 



where jo and ji are the spherical Bessel functions. 
Here we introduced the linear power-spectrum defined by 

(<5(ki)<5(k 2 )) = P(fci)<5r,(ki+k 2 ). Next we define at z = 0: 



vl = 4tt (H f(n mfi )y 



dkP{k;0) 



and 



Il v (x) = 3 



JdkP(k;0) \ ]Q (kx)-2^1 



JdkP{k;0) 



(12) 



(13) 



Thus, IL v (x) only depends on the shape of the power- 
spectrum and we have n„(0) = 1 while at large scales 
x 3> R—i it shows an oscillatory behaviour. Then, we can 
write £«( x ) as: 



Ux) = -£M*) 



0- + Z) 



_ 2 H(z) f[Q m (z)] 

Ho /(^m : o) 



(14) 



where we used the redshift evolution P(k) <X (1 + z) 2 of 
the linear power-spectrum. 



(8) 2.3. Free electron density fluctuations 



Second, we have to model the fluctuations of the free 
electron density distribution. As seen in (EJ) two effects 
contribute to these fluctuations: i) inhomogeneities of the 
baryonic matter distribution and ii) of the ionization frac- 
tion. In most previous studies of these secondary CMB 
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anisotropics only the second contribution was taken into 
account through a model of ionized bubbles within a uni- 
form IGM (e.g., Gruzinov & Hu 1998; Knox et al. 1998). 
However, as shown in Valageas & Silk (1999a) the clump- 
ing of the gas is not negligible even at z ~ 10. As can be 
seen from (0), and as we shall check below, this increases 
somewhat the amplitude of these CMB anisotropies. In 
this article we use a simple model to estimate the cor- 
relation term £ e (r). First, although the baryonic density 
fluctuations may be correlated with the ionization fraction 
we make the approximation: 

C e ~ ((1 + <5i)(l + 6 2 )){x e , lXe , 2 ) - (1 + S) 2 (x e } 2 (15) 

which can be written: 

Ur) = ([1 + 6GM(r)](Xe^e,2}W - (x e ) 2 ) (l+S) 2 GM .(16) 

Here we note £igm(0 the two-point correlation function 
of the gas density within the IGM and (S)igm the mean 
density contrast of the IGM. Note that at low redshifts 
after reionization x e = 1 almost everywhere within the 
IGM and the relation £ e (?") = £igm(?~)(1 + S) 2 GM becomes 
exact. Until z ~ 1 most of the matter is contained in the 
IGM since only a very small amount of gas was able to 
cool and form galaxies, hence (1 + <5)igm is close to unity, 
while at z — we have (1 + <5)igm ~ 0.4 (Valageas & Silk 
1999a; Valageas et al. 2000). 

2.3.1. Ionization fraction 

In order to obtain the term (x e , 12^,2) (?") we use a model of 
spherical ionized bubbles around galaxies and quasars as 
in Valageas & Silk (1999a). Thus, we consider that x e = 1 
within ionized patches and x e — everywhere else. At low 
z after overlap of the ionized regions x e — 1 throughout 
the IGM. Indeed, reionization occurs thanks to the growth 
of the ionized bubbles which finally occupy all the volume 
(and not through a slow increase of a uniform ionization 
fraction). Then, (x e ^x e ^){r) is simply the probability that 
two points at distance r are located within ionized regions. 
First, let us consider uncorrelatcd ionizing sources. Then, 
the probability that the first point ri is within an ionized 
bubble is the volume fraction occupied by these regions: 
da;i p 



x\ M\ 



(17) 



where Mon(^i) is the volume of the ionized bubble 
associated to the source x\ (galaxy or quasar) and 
(jj / Mi)xiH(xi)dx\ is the number density of radiation 
sources labelled by the parameter x\ in the range [x\, Xi + 
dxi] (see Valageas & Silk 1999a; Valageas & Schaeffer 
1997). Next, a simple geometrical calculation shows that 
the second point r 2 = ri + r can be located within the 
same ionized spherical bubble of radius R\ with a proba- 
bility V 2 ,s- 



V 2 



F 



Ri 



with 




ifx < 2 



On the other hand, the point r 2 can be embedded within 
another bubble with a probability 7^2, o with: 



(18) 



1-F 



Ri 



dx 2 p 
'x~ 2 ~W 2 



x 2 2 H(x 2 )V ion (x 2 ) (19) 



where we assumed that the bubbles associated to differ- 
ent sources do not overlap. This should be a good ap- 
proximation before reionization when there are very few 
ionized regions. Hence we obtain for uncorrelated sources 
(x et ix e , 2 )(r) = Vi(V 2jS + V 2 , )- However, the sources x\ 
and x 2 should be correlated, especially at high z where 
they correspond to very rare large density fluctuations. 
Thus, we write: 



(20) 



V 2 (r;xi)=Mm{l, F[ — 



1-F — 
Ri 



I —^-4H(x 2 )l\ on (x 2 ){l+^ UX2 (r)} 

x 2 M 2 



where S, Xl ,x 2 ( r ) is the correlation function of the objects 
x\ and x 2 . Following the analytic results obtained in 
Bcrnardeau & Schaeffer (1992, 1999) we have: 



£.x u x 2 (r) = b(xi)b(x 2 )t(r) 
with 



(21) 



i«l: b(x) ex x (1 ^ )/2 and x > 1 : b(x) ex x (22) 

where b(x) is the bias associated to an object defined by 
the parameter x. This behaviour has been shown to agree 
with the results of numerical simulations in Munshi et al. 
(1999). In Valageas et al. (2000) the predictions of this 
model for the correlation functions of galaxies, quasars, 
Lyman-a clouds and clusters have been compared to ob- 
servations. Thus, this provides a unified model for all these 
objects (which are characterized by different scales and 
densities) which agrees reasonably well with observations 
and simulations. In (E0[) the use of the minimum ensures 
that V 2 < 1. Indeed, at small redshift after reionization 
the second term in the minimum of (120) becomes larger 
than unity since all the volume is ionized (the ionized bub- 
bles overlap and Vi on goes to infinity). Finally, we obtain 
(xe,ix e , 2 )(r) as: 



(23) 



(x ei ix ei2 ){r) = 

Min/l , I ' ^ll-x\H{xi)V ion {xx)V 2 {r;xx) 

I J xi Mi 

while the mean ionization fraction is simply: 



,'■. = Mi"- < I ■ / ^-^-xlHix^Vionixi) } . (2-1 1 



The minima in (|23|) and (|24]) ensure again that 
(x e \X e 2 ) < 1 and (x e ) < 1 after reionization. 
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2.3.2. Density fluctuations 

Finally, we need to evaluate the two-point correlation 
function £igm(0 of the gas density within the f GM which 
appears in (|l6|). We consider that the universe is made 
of collapsed objects which have been able to cool and to 
form galaxies, embedded within a lower density medium 
which we call the IGM. Hence, the latter corresponds to 
voids as well as to density fluctuations (which may ap- 
pear as filaments or shallow spherical halos) associated 
with the Lyman-a forest (Valageas et al. 1999). Then, as 
in Valageas & Silk (1999a) the mean density of the matter 
which forms the IGM is given by: 



(1 + S) 1GM = (1 + A) u + / — x 2 H{x) 

Jo x 

while the mean square density is: 

r««oi dx 



(25) 



((1 + <5) 2 }igm = (1 + A)i + / — x z H(x){\ + A). (26) 

Jo x 

The term (1 + A) u corresponds to voids while the second 
term in the r.h.s. in ( |25|) and (|2q ) is the contribution from 
Lyman-a forest clouds. The upper bound x coo i ensures 
that we do not count cooled objects (which are identified 
to galactic halos) nor clusters of galaxies, since they are 
not part of the IGM. Then, the clumping factor of the 
IGM is simply: 



a 



IGM 



1 + £igmO" = 0) = 



((l + ^IGM 

(l + S) 



2 
IGM 



(27) 



(which we noted C„ in Valageas & Silk 1999a). In addition, 
as described in Valageas et al. (1999) and Valageas & Silk 
(1999a), baryonic density fluctuations within the IGM are 
erased on small scales below Rd(z) due to the non-zero 
temperature Tigm with: 



Rd(z) ~ t H \ 



' kT\GM 

fim p 



tH c s 



(28) 



This corresponds to the scale reached by acoustic waves 
over the Hubble time tn (with the sound velocity c s ). 
Note that the damping scale Rd is usually smaller than 
the standard Jeans scale -Rj ean s which corresponds to the 
limit of large times. Hence we write: 



r < Ra(z) : 6gm(0 = ^iGu(Rd) = Cigm - 1 



(29) 



in order to take into account the damping of the baryonic 
power-spectrum on small scales within the IGM, in a fash- 
ion which is consistent with (E7p. Finally, for larger scales 
we use the prescription: 



r > R d (z) : 6gmM = (Cigm - 1) ttttt 



(30) 



Note that this large— r behaviour (i.e. for r > Rd) is con- 
sistent with the factorization (Ell). 



2.4. Contribution from galactic halos 

In addition to the free electron IGM number density corre- 
lation function £ e (r) which we defined in (|16[), we also in- 
troduce the correlation function £ e , ga i if) due to collapsed 
halos which were able to cool and to form galaxies. We es- 
timate this contribution as follows. We assume that within 
these halos the gas is totally ionized (by the radiation of 
the central stars or QSO and by collisional ionization, due 
to shock- heating up to the virial temperature). Then, we 
can write (15l) as: 



£e, g ai(r) = F g 2 a , vol (1 + A c ) 2 £ gal (r) = F g 2 al £ gal (r) (31) 

where -Fgai.voi is the volume fraction (i.e. the filling fac- 
tor) occupied by these objects and F ga \ = F ga i )V0 i(l + A c ) 
is the fraction of matter they contain. We note A c (z) the 
mean density contrast of just-virialized halos, given by the 
usual spherical model, while £ ga i(r) is the matter correla- 
tion function associated with these halos. If we do not 
take into account the density profile of these objects (i.e. 
we take a constant density equal to (1 + A c )p) we have for 
the clumping of the gas: 



C g al = l+£gal(r = 0) = 



1 



F, 



>1. 



gal 



(32) 



Then, in a fashion similar to (|29j) and ( |3C| ) we write the 
matter correlation function as: 

£(r) 



r > Rcoo\(z) : £gal( r ) = (Cgal - 1) 



(33) 



eoRcooO 

and £ g alO") = (Cgai — 1) for r < R coo \- Here we defined 
^cooi(z) as the virial radius of the smallest halos which are 
able to cool at redshift z, as in Valageas & Silk (1999a). 
Of course, our prescription for £ e ,ga,\(r) is only meant to 
provide a crude estimate of the contribution due to galac- 
tic halos. Note that if we consider a specific density profile 
for the dark matter halos we would only need to multiply 
Cgai by a numerical factor of order unity, provided that 
the slope of the density in the inner regions is shallower 
than —3/2 (otherwise C ga i is not finite). In addition, one 
should take into account the additional collapse of the gas 
(while it cools) and the effects of matter ejection through 
supernovae or stellar winds. In this article, we shall restrict 
ourselves to the simple prescription ( p3| ) which should ex- 
hibit the right trend with r, keeping in mind that at small 
scales r < -R C ooi we might underestimate the signal by a 
factor ~ 3. The comparison of £ e . ga i with £ e allows us to see 
which of these two contributions (from galactic halos or 
from the IGM) is dominant at a given scale. It also shows 
the characteristic scales associated with both processes. 

2.5. Power spectrum, of secondary anisotropies 

Using the results of the previous sections, we can write the 
two-point correlation function C(9) obtained in un) as: 

dr rn„[r(l + z)R e (r) 



C(6) 



2tX 
3 c 2 



/ dz w(z) / — 

J jRg K<Z 



V^— R2 



-(34) 
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with: 



w(z) 



and: 



dX 



V(z) 



H(z)f[n m (z)] 



c V(z) 
2 H 1 + z' 



(35) 



(36) 



In (|34|) we made the change of variable r = |Rg + 1|. The 
angular power-spectrum of these secondary anisotropies is 
given by: 



Ci=2ir d(cos(9) Pi(cos9) C{6) 

which for I 3> 1 can be approximated by: 

I > 1 : Q = 2tt / d9 J (W) C(6) 



(37) 



(38) 



where P; are the Legendre polynomials and Jo is the Bessel 
function of order 0. Then, from (B4j) we obtain after inte- 
gration over 9 (see Gradshteyn & Ryzhik 1965, §6.554.2, 
p.706): 



C, 



^ 2 v 2 
3 c 2 I 



dz w(z) I -£- -£- sin f -£- 



Rz Rz 



Rz 



xn„[r(l + Z )]^ e (r). (39) 

The expression (p9) clearly shows that C\ will exhibit 
a strong decline at large I for wavelengths which are 
smaller than the typical scales of the free electron in- 
homogeneities (R z /l <C Rn c ) because of the oscillatory 
term sm(l.r/R z ). On the other hand, for small I (i.e. at 
large scales R z /l ^> R Ua where the correlation of the elec- 
tron distribution vanishes) we recover a white noise power- 
spectrum (Ci is constant). 

Finally, following Bruscoli et al. (2000) we define the 
quantity Si by: 



* = i 



31/2 



dllCi 



1/2 



, /(/ + l)Ci 

* 2n 



(40) 



This "local average" of 1(1 + 1)Ci/(2tt) removes some of 
the oscillations of Ci at large I and it allows us to see more 
clearly the drop at large / of the spectrum. It is also more 
convenient for numerical calculations since it requires a 
lower resolution. From the relation (B9|) we obtain: 



Si 

Using 
S,= 



2t X 



3 c 2 



dz w(z) 



dr 

R~z 



Ir 
2R Z 



Sir 
2R Z 



xiyr(l + z)K e (r). (41) 
we can also write Si in terms of C{9) as: 



de 



C(9) 



316 



-Ji 



3W 



19 



-•/: 



(42) 



Note that the expressions (|I|) , ( J39J ) and (^Tj) are quite 
general and do not depend on our model for the inho- 
mogeneities of the free electron distribution (which enters 
£ e (?"))• They merely use the linear evolution (Q) of the 
velocity fluctuations and the approximation that velocity 
and density fluctuations are uncorrelated. Moreover, it is 
more accurate (hence more convenient in terms of numer- 
ical resolution) to compute Ci and Si from (|39"|) and (41) 



rather than from the Fourier transform of the final angular 
correlation C(9). 

3. Numerical results 

For the numerical calculations we consider an open CDM 
universe (OCDM) with fl m = 0.3, Ov = 0, n h = 0.03, 
H n = 60 km/s/Mpc and as — 0.77. These values are 
those we used in previous articles where we considered 
the luminosity functions of galaxies (Valageas & Schaeffer 
1999), Lyman-a absorbers (Valageas et al. 1999), clusters 
(Valageas & Schaeffer 2000) and reionization by stars and 
quasars (Valageas & Silk 1999a, b). 

3.1. Velocity fluctuations 

First, we display in Fig.pl the magnitude of the velocity 
correlation function £„(#) at z — 0. It shows a plateau at 
small scales x <§C R-i (most of the power comes from i?_i 
where the quantity kP(k) is maximum) and it declines at 
large scales with oscillations. Thus, in our model velocity 
fluctations are of order Av ~ 400 km/s at z = (and 
they roughly decrease as (1 + z) -1 / 2 ). They only arise 
from the linear growth of initial perturbations through 
gravitational interaction. In the actual universe, an ad- 
ditional source of velocity fluctuations (mainly at small 
scales) could be provided by other processes like the ejec- 
tion of matter by supernovae or turbulence. If these ve- 
locity fluctuations are larger by a factor a than the value 
we use for Aw shown in FigjU then on the correspond- 
ing scales we should roughly increase our predictions for 
the secondary CMB anistropies Ci by a factor (1 + a 2 ). 
We shall not investigate here this possibility but we can 
note that a significant effect would require rather large 
velocities. 

3.2. Free electron correlation junction 

The angular correlation C{9) and the power-spectrum Ci 
correspond to the sum over the line of sight of the fluc- 
tuations of the free electron number density, up to the 
recombination redshift. Then, this integration over red- 
shift could blur some features of these density fluctuations. 
Hence, in order to clarify the analysis it is interesting to 
consider first the real-space correlation function £ e (r;z) 
obtained for a given redshift z. This also allows us to see 
the evolution with redshift of the free electron density fluc- 
tuations. 
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Fig. 1. The magnitude of the velocity correlation function 
£, v (x) at redshift z — 0. The solid curve shows the quan- 
tity v \a v (x)\ 1 / 2 , see Q. The feature at x ~ 100 Mpc 
corresponds to the fact that H v (x) becomes negative for 
X > 100 Mpc. 

In order to understand the physical origin of the sig- 
nal we split up the correlation function into several parts. 
First, we consider the total contribution J e (r) from the 
IGM, as defined in ([lq). Next, we introduce £ £]U as the 
correlation function we get when we do not take into ac- 
count the correlations of ionized bubbles (subscript "u" for 
"uncorrelated" ) . That is, in ( J20| ) we set the term £,x 1 ,x 2 { r ) 
to 0. Then, we define £ e ,h as the signal produced by un- 
correlated bubbles into a "homogeneous" IGM (subscript 
"h" for "homogeneous"). That is, in (n6|) we set the term 
(igmW to 0. Thus, £ ej h allows us to see the contribution to 
£ e due to the inhomogeneities of the free electron number 
density produced by patchy reionization in distant bub- 
bles. Then, £ ejU shows by comparison with £ ej h the impor- 
tance of the clumping of the gas within the IGM. Finally, 
the difference between £ ejU and the total signal ^(r) mea- 
sures the effect of the correlations of these ionized bubbles. 

We show our results in Fig.|| for two different red- 
shifts. The upper panel at Z\ = 7.5 corresponds to a 
redshift slightly before reionization (at z ri = 6.8). First, 
we clearly see on £ ei h the contribution due to patchy 
reionization within finite size bubbles. Thus, £ e .h is con- 
stant at small scales below the characteristic size of the 
ionized bubbles (~ 0.5 comoving Mpc) and it drops at 
large scales. Of course, since £ e ,h only corresponds to fluc- 
tuations in the ionized fraction (it does not take into 
account matter density fluctuations) it is of the form 
£e,h(»") = (xe,i x e,2)(f) — (x e ) 2 . Since by definition we have 
x e < 1 we get £ ei h < 1- In fact, £ ej h never reaches unity be- 
cause at large redshift the filling factor Qhii of the ionized 
bubbles is much smaller than unity while at low z after 
reionization all the medium is reionized hence £ e ,h CO = 0. 
In the upper panel we have Qhii — 0.9 < 1 while in the 
lower panel after reionization we have £, e .h( r ) = hence 
the curve does not appear in the figure. 
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Fig. 2. The real-space two-point correlation function £ e (r) 
at two different redshifts. The solid line labeled £ e (resp. 
£e,gai) shows the contribution from the IGM (resp. from 
galactic halos). The curves £ ejU and £ ej h correspond to the 
"uncorrelated bubbles" and "homogeneous IGM" scenar- 
ios (see main text). 

Next, £ 6]U shows the influence of the clumping of the 
gas within the IGM. Since Cigm > 1 we no longer have the 
upper bound £ eiU < 1 and we can check in the figure that 
indeed we can have £ ejU > 1. However, the typical overden- 
sities within the IGM are smaller than the density contrast 
A c (z) of just-virialized halos hence Cigm < 1 + A c , see 
the expression (^7|). Thus, we obtain £ EjU < A c (z). We can 
check in both panels in Fig.|] that our results agree with 
this upper bound (A c (z) ~ 200). In fact, £ EjU is signifi- 
cantly smaller (£ 6]U ^ 10) since most of the volume of the 
universe is filled by lower density regions with p <, ~p. In 
agreement with (Ey) the correlation function £, e ,u{f) satu- 
rates at small scales below the damping scale i?d and it fol- 
lows the decrease at larger scales of the matter correlation 
function. Since Rd is smaller than the typical sizes of the 
ionized bubbles the characteristic break of £ e! u(f) occurs 
at smaller scales than for £ e! h- In other words, the mat- 
ter density fluctuations provide more small scale power in 
relative terms than ionized bubbles. Moreover, since the 
clumping of the gas is rather large, even at large redshifts 
(at z > 10 we already have Cigm ^10, see Valageas & 
Silk 1999a), we find that £ ejU » £ ej h at all scales. 

However, the presence of the ionized bubbles is not 
totally blurred by the superimposed matter density fluc- 
tuations. Indeed, we can see in the upper panel that at 
large scales the actual correlation £ e is much larger than 
£e !U . This means that for r ^ 1 comoving Mpc the signal 
is dominated by the cross-correlation of ionized bubbles. 
This arises from the correlations of their central collapsed 
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halos (which are the sites of formation of the central ra- 
diation source, either a galaxy or a QSO), see (p0|), which 
correspond to very rare overdensities. Of course, by defi- 
nition this effect appears at scales larger than the typical 
size of the ionized bubbles. It provides excess large scale 
power to £ e , as compared with £ e>u , and it leads to a char- 
acteristic feature in the shape of the correlation function 
£ e (r). After reionization (lower panel) since there are no 
more ionized bubbles this effect disappears and £ e (r) is 
only governed by the fluctuations of the matter density, 
hence £ ejU (not shown) becomes equal to £ e - 

Finally, the curve £ e ,gai shows the contribution from 
galactic halos. We can see that the scale i? coo i is some- 
what smaller than the damping scale i?d (this is related 
to the fact that the density of virialized halos is larger by 
a factor A c ~ 200 than the typical IGM density) so that 
£e,gai saturates at smaller scales than £ e . Moreover, the 
clumping factor associated with these halos is larger than 
for the IGM (because of this difference between the typ- 
ical densities) hence the plateau at small scales of £ e , ga i 
is higher than for £ e . From (p^) and ((32|) we can check 
that £e,gai < F ga iA c < A c . We can see in the figure that 
£e,gai is somewhat smaller than this upper bound because 
the fraction of matter enclosed within such halos is still 
low at these redshifts: F ga i ^ 0.12 at z^ 7 (see also 
Fig. 12 in Valageas & Silk 1999a). We note that at larger 
scales (r ~ 1 Mpc) the contribution from the IGM and 
"galaxies" are of the same order while at very large scales 
(r > 1 Mpc) before reionization the signal is dominated by 
the IGM through the cross-correlation of ionized bubbles. 
We can expect to recover these features in the integrated 
quantities C(6) and Ci. 
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Fig. 3. The angular two-point correlation function C(8). 
The solid line labeled C shows the contribution from the 
IGM. The curve C ga i displays the contribution from galac- 
tic halos. The curves C u and Ch correspond to the "uncor- 
rected bubbles" and "homogeneous IGM" scenarios (see 
main text). 

and C over all angles 8 (in particular down to 8 — > 0) 
since a small angular separation 8 corresponds to a large 
real-space distance r at high z. This leads to a difference 
with the real-space correlation £ e (r) shown in Fig.0 where 
we found that at small scales (below the size of ionized 
bubbles) £ e-u is equal to £ e - In a similar fashion, the inte- 
gration over redshift also leads to smoother curves C(8). 
Finally, we can see that the contribution C ga i from galac- 
tic halos is larger than the signal from the IGM at small 
scales 8 Ss 10~ 4 rad while it becomes smaller at larger 
scales, as expected from Fig.g. However, on the whole the 
difference between both contributions is not very large. 



3.3. Angular two-point correlation function 

We show in Fig.H our results for the angular two-point cor- 
relation function C(8). We can check that we recover the 
trends described in Sect. 3. 2, which have not been totally 



blurred by the integration over redshift. 

First, we recover as in Fig.g the characteristic shape of 
the contribution Ch from pure ionized bubbles: a plateau 
at small angular scales 8 <2 10 -4 rad and a sharp drop 
at larger scales, beyond the size of the bubbles. Note that 
a characteristic length scale r is related to the angular 
scale 8 by r ~ 8c/ Ho ~ 58 x 10 3 (in comoving Mpc). 
Hence 8 ^ 10~ 4 corresponds indeed to the scale r ~ 0.5 
Mpc seen in Fig.0 (the knee of £e,h(f))' Then, the correla- 
tion C u shows a smoother shape, due to the power-law 
behaviour of the real-space matter correlation function 
£(r), with a break at a smaller scale 8 ~ 10~ 55 rad due 
to small scale matter density fluctuations. Moreover, we 
have C u > Ch, in agreement with Fig.H Next, the total 
signal C(8) from the IGM is larger than C u , especially at 
large angular scales 8 ~ 10~ 3 rad, because of the cross- 
correlation of ionized bubbles. Note that the integration 
along the line of sight spreads the difference between C u 



3.4- Power- spectra Ci and Si 

We show in Fig .§ the quantity l(l + l)C[/(2n) and the "lo- 
cal average" Si, for the contributions from the IGM and 
from galactic halos. We also plot for comparison the power 
spectrum of primary anisotropies calculated using CMB- 
FAST (Seljak & Zaldarriaga 1996). At small I both quan- 
tities 1(1 + l)Ci/(2ir) and Si are almost identical since Ci 
varies slowly with the wavenumber I. At large I (I ^ 10 5 ) 
the power-spectrum Ci exhibits an oscillatory behaviour 
(since it is the Fourier transform of a function C(ff) which 
shows a sharp drop at large angular scales) and a slow 
decline. In particular, at very large I (I ^ 10 7 ) the oscilla- 
tions of Ci are not resolved by the numerical calculation. 
On the other hand, Si becomes significantly different from 
1(1 + l)Ci/(2ir) as it shows a sharp decrease with I and 
fewer oscillations. Of course, this is due to the "averaging 
procedure" which enters the definition (BTJ) of Si. At large 
I the numerous oscillations of Ci over the range [1/2, 31/2} 
partially cancel out which leads to a stronger falloff for 
Si. Moreover, this "averaging" smoothes the behaviour of 
Si which shows much weaker oscillations. This also allows 
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us to resolve Si up to larger I (note that Si is not com- 
puted from Ci but directly from the expression (El)). This 
suggests that for observational purposes too, the quantity 
Si may be more convenient as it should be more robust 
(i.e. require a lower resolution) than Ci at large I and it 
shows more clearly the transition to the large-/ regime by 
a sudden drop. 

Note that for a correlation function (7(0) = 
exp[— (0/0o) 2 ] with a Gaussian cutoff with a characteristic 
scale 0q we get from (E3): 



-(W Y/1G 



while (pq) gives: 



2tt 



(Wo) 



-9(;e f/i6 



9o) 2 /4 



e f/i6 



(43) 



(44) 



On the other hand, for an exponential cutoff (7(0) = 
exp(— 0/0o) we obtain the same slow cutoff for both quan- 
tities: Si ~ l 2 Ci ~ 1/(19q). Finally, for a correlation func- 
tion which is a top-hat (i.e. (7(0) = 1 for 9 < 9q and 
(7(0) = for 9 > So) we get: 



Si = Jo 



W 



-Jo 



319q 



and 



2tt 



(W )J 1 (Wo) 



(45) 



(46) 



In this latter case, at large I both (7/ and Si show an 
oscillatory behaviour but we have \Si\ ~ Z -1 / 2 while 
\l 2 Ci\ ~ Z 1 / 2 , so that Si shows again a stronger decrease 
than 1(1 + 1)Ci/(2tt). Thus, these examples explicitly show 
that the shape of the large-/ tail of Si depends rather 
strongly on the details of the angular correlation function 
(and even more so for C{). However, the "local average" Si 
usually shows a drop beyond the characteristic wavenum- 
ber l/0o (while \l 2 Ci\ may either decrease or grow in an 
oscillatory fashion). 

The oscillatory behaviour of Ci (hence the sharp drop 
of Si) appears at lower I for the IGM than for galactic 
halos. This is due to the fact that the correlations £ e {r) 
and (7(0) show less small-scale power for the IGM con- 
tribution than for the signal from galactic halos, as seen 
in Fig.|| and Fig.||. Indeed, as shown in (38) and ( |39| ) a 



characteristic angle (resp. a physical length r) translates 
into a characteristic wavenumber I ~ 1/0 (resp. I ~ R z /r). 
Hence, since the scales of the fluctuations of the free elec- 
tron number density are larger for the IGM than for the 
galactic halo component the power-spectra Ci and Si of 
the IGM appear shifted towards smaller I with respect to 
the contribution from galaxies. At low I we recover a white 
noise spectrum (Ci is constant, hence Si oc I 2 ) since this 
corresponds to very large scales where the correlations of 
the electron distribution are negligible. However, the slope 
of Si we find at low I for the contribution from galactic 
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Fig. 4. The power-spectra Ci and Si of the secondary 
anisotropies for the OCDM cosmology. The solid curves 
show the quantity Si for the contributions from the IGM 
and from galactic halos. The dotted (resp. dashed) curve 
with oscillations at large I displays I (I + 1)Ci/(2tt) for the 
IGM (resp. galactic halos). The upper curve labeled "pri- 
mary" shows I (I + 1)Ci/(2tt) of the primary anisotropies 
for the same OCDM model (see text). 

halos is smaller because of the large-scale correlations of 
these rare overdense objects. 

We note that we clearly recover the main features of 
the correlation (7(0) shown in Fig.g. Thus, for the IGM 
we find that the transition to the white noise part (the 
falloff at low I with a slope I 2 ) occurs at / ~ 10 3 which 
corresponds to the cutoff at large angles ~ 10 -3 rad of 
the correlation (7(0) (see the strong knee in Fig.0). On the 
other hand, for the galactic halo contribution the transi- 
tion to the large-/ regime (marked by the drop of Si) ap- 
pears at I ~ 10 6 which corresponds to the scale ~ 10~ 6 
rad below which (7 ga i(0) saturates (hence to the smallest 
angular scale of the density fluctuations). Note that addi- 
tional power at smaller scales (due to substructures within 
halos and to the subsequent collapse of baryons when they 
cool) would shift this transition towards higher /. 

In a fashion similar to Fig.0 we can split up the spec- 
trum Si into several components. This decomposition is 
shown in Fig.pl First, in agreement with Fig.0 and Fig|3 
we can see that the contribution Si^ from uncorrelated 
ionized bubbles within a uniform medium is strongly 
peaked at I ~ 10 4 which corresponds to the typical size 
of the ionized bubbles (0o ~ 10 -4 rad, r ~ 0.5 comov- 
ing Mpc). At small wavenumber (Z0o) <^ 1 we recover a 
white noise spectrum Si ~ (/0o) 2 while at larger I we get 
a somewhat smoother decrease than l~ 2 . 

Next, in agreement with Fig.|| we can see that Si^ u is 
larger than S 1 /^ and it shows a broader maximum, due to 
the additional power at small and large scales provided 
by the density fluctuations within the IGM. The com- 
parison with the total signal Si shows that this Fourier 
transform actually separates the various contributions to 
the power-spectrum. Thus, for I > 10 4 we have Si = Si u 
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Fig. 5. The power-spectrum Si of the secondary 
anisotropies for the OCDM cosmology. The solid curve 
labeled Si (resp. 5/ )ga i) shows the contribution from the 
IGM (resp. from galactic halos). The curves S; iU and S^h 
correspond to the "uncorrected bubbles" and "homoge- 
neous IGM" scenarios. 

and the peak at I ~ 10 55 corresponds to the break at 
9 ~ 10 -5 " 5 of the correlation C u (8) (and to the saturation 
at r ~ 2 x 10~ 2 Mpc of £ e , u (r)). This directly probes the 
small-scale fluctuations of the matter density field. More- 
over, the equality Si = Si, u in this large-Z regime translates 
the fact that at small scales r < 1 Mpc we had £ e = £ eu , 
as seen in Fig.[| Note that in the angular space represen- 
tation C(9) the integration over redshift along the line of 
sight destroys this feature as C{9) > C n {9) for all angles 
and one cannot recognize from the total signal C the sig- 
nature of this small scale feature. Thus, the Fourier trans- 
form presents the strong advantage to separate various 
physical processes as they act on different scales. However, 
in this regime I > 10 4 the signal should be dominated by 
the contribution from galactic halos. 

Here we can note that in our calculations we modelled 
ionized regions as spherical bubbles while detailed numer- 
ical simulations show they can display a more complex 
morphology (e.g., Abel et al. 1999) as ionization fronts 
propagate more easily in voids. This means that the bell- 
shaped curve Si t h in Fig|| (which measures the contri- 
bution from patchy reionization, i.e. from the geometry 
of HII regions) underestimates the actual signal at large 
wavenumbers (I > 10 ) where we neglected the contribu- 
tion from higher-order spherical harmonics (which provide 
some power on scales smaller than the typical radius Ri 
of the ionized bubble). However, we can check in Fig.0 
that even if we spread the maximum of the curve S^h up 
to wavenumbers / ten times larger (i.e. the geometry of 
the ionization front displays significant power up to scales 
ten times smaller than Ri) our results remain unchanged. 
Indeed, most of the power in this range is due to the fluc- 
tuations of the matter density itself rather than to the 
geometry of HII regions. Moreover, the smallest scales dis- 
played by the geometry of the ionization fronts are at least 



of the same order as the size of the smallest structures of 
the density field (from which they originate). Then, the 
factor ri 2 in the correlation C(9) ensures that if a large 
fraction of such clouds, enclosed within the distance Ri to 
the radiation source, are ionized and if they show a density 
contrast 8 i<S 10, the signal is dominated on these scales by 
the fluctuations of the free electron density rather than by 
sheer geometrical patterns. Hence our results are not very 
sensitive to the approximation of spherical HII bubbles. 

Then, the comparison of Si with Sj iU shows that the 
peak at I ~ 10 3 comes from the correlations of ionized 
bubbles. This also corresponds to the knee at 6 ~ 10 -3 
rad of C{9). In agreement with Fig. g we find that at small 
wavenumbers I ^ 10 3 the power-spectrum is dominated 
by the IGM contribution (i.e. the correlations of ionized 
bubbles) while at larger I most of the signal comes from 
galactic halos which provide much small scale power. The 
peak at I ~ 10 6 of Sj )ga i corresponds to the break at 9 ~ 
10~ 6 rad of C ga \. Thus, we see that the power-spectrum 
of the final signal mainly probes the small-scale density 
fluctuations which give rise to galactic halos (and their 
possible substructures) and the correlations at large scales 
of ionized bubbles. 

Analytical studies of the power-spectrum Ci of these 
secondary CMB anisotropies have also been presented in 
some earlier works. Thus, Gruzinov & Hu (1998) consider 
the effect of patchy reionization within a uniform medium, 
assuming a Gaussian falloff for the real-space correlation 
function £ e (r). This could arise from a Gaussian distribu- 
tion of the size of the HII regions. Then, as in (E3) and 
(E3) they obtain a Gaussian cutoff for the angular corre- 
lation C{9) and the power-spectra Si and Ci. This corre- 
sponds to our curves Ch and Si t h in Fig.|| and in Fig.|[ 
Our predictions for this scenario are similar to their re- 
sults, but our spectrum Si peaks at a larger / (/ ~ 10 4 ) 
than theirs (I ~ 10 3 ). This is due to the fact that in our 
model the comoving size of HII regions is of order 0.5 
Mpc (see Fig.||) while they assume a very large radius of 
20 Mpc for the bubbles. We can also note that our cut- 
off at large I is smoother than a Gaussian. Indeed, since 
observed luminosity functions usually show a simple expo- 
nential cutoff (and our results match observations at low 
z) we can expect a shallow cutoff of the form exp(— r 1 ' 3 ) 
(because the volume of ionized bubbles is proportional to 
the number of ionized atoms, hence to the luminosity of 
the source) and we noticed above that a pure exponential 
cutoff already leads to a simple power-law decline of the 
spectrum Si (as I /I). Using a slightly more sophisticated 
model, Aghanim et al. (1996) calculate the reionization 
from early formed quasars, deducing the statistic of the 
ionized bubbles size from the distribution of quasar lumi- 
nosities. Their results predict that most of the power is at 
I ~ 10 3 and they are similar to a one-patch scenario with 
a bubble radius of ~ 10 Mpc, except for the high I cutoff 
which is less steep, due to the smaller patches distribution. 
Then, Knox et al. (1998) consider the effect of the corre- 
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lations of these ionized bubbles. In agreement with our 
results, they find that this leads to a much broader distri- 
bution of the power-spectrum. Note that in our analysis 
we have split up the influence of matter correlations into 
two processes: the cross-correlation of HII regions them- 
selves (through the correlation of the emitting sources), 
which provides additional power at larger scales (I < I0 4 ) 
than the size of these patches, and the fluctuations of the 
matter density field within these bubbles, which builds 
power at smaller scales (I > I0 4 ). Thus, the final signals 
C(9) and Si are the sum of the contributions from vari- 
ous processes. The advantage of our approach is that it 
provides a fully consistent description of these different 
scales, from a model built to study the detailed reioniza- 
tion history of the universe which has already been com- 
pared with observations for various aspects (e.g., galaxy 
luminosity function in Valageas & Schaeffer 1999; X-ray 
emission from clusters, galaxies and quasars in Valageas 
& Schaeffer 2000). 

3.5. Redshift distribution 

The angular correlation C(9) and the power-spectra Ci 
and Si correspond to an integration along the line of sight 
of the fluctuations of the free electron number density. 
However, it would be interesting to see the relative im- 
portance of the contributions from various redshifts to the 
final signal. In particular, this would show whether these 
secondary CMB anisotropies arise from a narrow range 
of redshifts close to reionization at z,i = 6.8 or from a 
more extended interval. Thus, we define the normalized 
quantity Cg(z) by: 



C e (z) 



^v 2 
ic 2 C{9) 



«,(*) / £ rlMr(l + z)U e (r) (4?) 



R z 



xA 



Rn 



It obeys the normalization condition: 

dzC {z) = 1 



(48) 



and Cg(z)dz is the fraction of the final angular correlation 
function which is generated in the redshift interval Az. In 
a similar fashion we define the quantities Ci{z) and Si{z), 
from (p9|) and pi]), which are also normalized to unity. For 
these three redshift distributions we alternatively consider 
the contributions from the IGM and from galactic halos. 

3.5.1. IGM contribution 

First, we consider the contributions from the IGM to the 
secondary CMB anisotropies. We show our results for the 
redshift distribution Cg(z) of the angular correlation in 
Fig.@, for two different angular scales. We note that the 
contributions to the final signal C(9) come from a rather 
large range of redshifts, typically 7.5 < z < 10 (so that 
Sz/Zxi ~ 0.36). There is a sharp cutoff at z ~ z r ; since at 
lower redshifts there are no more ionized bubbles. This 
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Fig. 6. The redshift distribution Cg(z) of the angular cor- 
relation (normalized to unity) from the IGM, for 9 = I0~ 3 
rad (solid line) and 9 = I0~ 5 rad (dashed line). 
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Fig. 7. The redshift distribution Ci(z) of the power- 
spectrum Ci (normalized to unity) from the IGM, for 
/ = 10 3 (solid line) and I — I0 5 (dashed line). 

drop is sharper for larger angular scales, in agreement 
with Fig.|| and Fig.^ where we noticed that large scales 
9 > I0 -4 rad (I < 10 4 ) are dominated by the correlations 
of ionized bubbles. However, at small scales 9 ^ 10~ 5 there 
is a non-negligible tail at lower redshifts due to the mat- 
ter density fluctuations within the fully ionized IGM. Of 
course, smaller scales also show a slightly more extended 
tail at high z since at higher redshift the typical size of 
ionized bubbles and the correlation length of the matter 
density field were smaller, which damps the contribution 
to large angular scales 9. Hence the redshift distribution 
of the angular correlation G(9) is somewhat broader for 
lower 9 (which translates into the smaller height of the 
maximum of Cg(z) in the figure since the curves are nor- 
malized to unity). 

We display in Fig.FI] the redshift distribution Ci(z) of 
the contribution from the IGM to the power-spectrum Ci 
for two different wavenumbers (normalized to unity). Of 
course, for I = 10 3 we recover a shape similar to the red- 
shift distribution Cg(z) we obtained for 9 = 10~ 3 rad, 
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Fig. 8. The redshift distribution Si(z) of the power- 
spectrum Si (normalized to unity) from the IGM, for 
I = 10 3 (solid line) and I — 10 5 (dashed line). 

since both quantities correspond to the same scale. For 
I = 10 5 the envelope of Ci agrees again with the shape 
we got in Fig.^J for C${z) (with a tail at low z due to 
small-scale density fluctuations) but the distribution Ci{z) 
now shows several oscillations. This is due to the Fourier 
transform involved in the definition of Cj. Thus, the con- 
tributions from successive redshifts along the line of sight 
almost cancel out. This agrees with the behaviour we ob- 
tained in Fig.0. Note that the oscillations occur before 
reionization: they are due to the patchy pattern of reion- 
ization in HII bubbles with a size larger than R z /l. At 
lower redshift this feature disappears as there are no more 
ionized regions to single out a large characteristic scale. 

Finally, we display in Fig.0 the redshift distribution 
Si{z) of the contribution from the IGM to the power- 



spectrum Si. As expected, for I 



10 3 



we recover 



the 



results we obtained for the power-spectrum C/. Indeed, 
as noticed in Figjj for small wavenumbers there are no 
oscillations since one probes scales which are larger or of 
the order of the correlation lengths of the free electron 
distribution, so that Si ~ l(l + l)Ci/(2n). At larger I some 
oscillations start to appear and Si(z) shows a different 
shape than C;(z). In particular, the oscillations of Si(z) 
are much smoother and broader than for Ci(z) and they 
appear at a larger wavenumber. Indeed, the "local aver- 
aging" over I associated with the procedure used to define 
Si "smoothes" the contributions from various scales. In 
particular, this allows us to see more clearly the redshift 
distribution associated with I = 10 5 where there are no 
oscillations yet. Moreover, it clearly shows that the large 
oscillations we obtained shortly before reionization for C'i 
almost cancel out so that high redshifts z > z T \ only pro- 
vide a small contribution to the final signal. This leads 
to a redshift distribution which is very different from the 
one obtained for smaller I which shows a sharp cutoff at 
z ~ z x \. Thus, we find that for these small scales the contri- 
butions to the power-spectrum Si come from an extended 
range of redshifts 2 < z < 8. As noticed in Fig.g, we find 



again that the use of the power-spectrum allows one to 
clearly see the various processes associated with different 
scales, which are somewhat blurred in the angular repre- 
sentation C(8). 

3.5.2. Contribution from galactic halos 
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Fig. 9. The redshift distribution C$ tga \(z) of the angular 
correlation (normalized to unity) from galactic halos, for 
6 = 1CT 3 rad (solid line) and = 10~ 5 rad (dashed line). 
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Fig. 10. The redshift distribution Si tga ,\(z) of the power- 
spectrum Si (normalized to unity) from galactic halos, for 
I = 10 3 (solid line) and I = 10 5 (dashed line). 

Now, we consider the redshift distribution of the con- 
tribution from galactic halos to the CMB anisotropies. We 
show in Fig.S our results for the angular correlation func- 
tion. First, we note that there is no drop at the reioniza- 
tion redshift z T \. Indeed, the fact that ionized bubbles sud- 
denly overlap so that the signal from the patchy pattern 
of reionization disappears does not affect the contribution 
from galactic halos. On the other hand, reionization does 
not lead to a sharp drop of the galaxy or quasar multi- 
plicity functions either since it does not imply a sudden 
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increase of the IGM temperature and of the Jeans mass. 
Indeed, as seen in Valageas & Silk (1999a) most of the 
reheating of the universe occured earlier in a gradual fash- 
ion so that the small increase of the IGM temperature at 
z T i has no impact on the population of radiation sources. 
Thus, the redshift distribution Cg iga x(z) follows the growth 
of non-linear structures so that smaller redshifts provide 
a larger contribution. This appears clearly for 9 = 1CP 3 
rad where most of the signal is generated at z ^ 1 when 
the scale r ~ 5 Mpc enters the non-linear regime. On the 
other hand, for the smaller angular scale 9 = 10~ 5 rad the 
contribution from very low z becomes smaller as the typi- 
cal size of virialized objects becomes larger than the scale 
which corresponds to 0. However, we may underestimate 
the power at low z because we neglected substructures 
within halos. Note that on these small angular scales the 
secondary CMB anisotropies should be dominated by the 
contribution from galactic halos. Hence they arise from a 
very broad range of redshifts (typically < z < 7) which 
is not related to z T {. 

We display in Fig.[l(] the redshift distribution Si >sa i(z) 
of the power-spectrum Si tga ,\- We recover a behaviour simi- 
lar to Fig|| In particular, note the large range of redshifts 
which is probed by large wavenumbers I ~ 10 5 . Higher I 
which are beyond the cutoff of the spectrum Si show in- 
creasingly important oscillations. 

3.6. Quasars versus stars 

In usual scenarios the universe is reionized by the radiation 
emitted by non-linear structures as collisional ionization 
is likely to be less efficient (e.g., Madau 2000; Valageas 
& Silk 1999b). There are two natural sources of radiation 
in present cosmological models: stars and quasars. In par- 
ticular, in our model the universe is reionized when HII 
bubbles created by galaxies and quasars overlap at z r i (see 
Valageas & Silk 1999a). The multiplicity functions we use 
for galaxies and QSOs are normalized to the low-redshift 
universe (z < 4) and are obtained in a consistent fashion 
(see also Valageas & Schaeffer 1999). Then, we find that 
the energy output provided by QSOs is of the same or- 
der as the energy radiated by stars. However, the spatial 
features of these two reionization processes may be differ- 
ent since one can expect QSOs to create fewer but more 
extended HII bubbles, since quasars are not as numerous 
as galaxies but their luminosity is much larger. Hence, the 
correlation function C{9) and the power-spectrum Ci may 
show more large-scale power for a quasar-driven reioniza- 
tion than for a galaxy-driven process. This is quite in- 
teresting as it might allow one to discriminate both sce- 
narios. Note on the other hand that the reionization of 
helium is usually due to the radiation emitted by quasars, 
as in our model. Indeed, stars have a black-body spectrum 
which yields very few high energy photons while quasars 
exhibit a harder power-law spectrum over the relevant fre- 
quency range. However, as pointed out by Tumlinson & 



Shull (2000) population III metal-free stars have a harder 
spectrum than typical low-z stars, hence they might be 
able to ionize helium in addition to hydrogen. Thus, the 
relative importance of quasars and stars is still an open 
problem. Unfortunately, we shall see below that the obser- 
vation of the secondary anisotropies of the CMB is unlikely 
to answer this problem. 
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Fig. 11. The angular two-point correlation functions 
C s (6) and Cq(8) due to stars and quasars. The solid line 
labeled C s (resp. Cq) shows the contribution from the 
IGM when we only take into account the ionized bubbles 
created by a central QSO (resp. by a central galaxy). The 
curves Ch, s (resp. C^q) corresponds to the "homogeneous 
IGM" scenario of ionized bubbles created by stellar radia- 
tion (resp. quasars) within a uniform medium. The dotted 
line shows the total IGM correlation function C(9) (it is 
identical to the solid line C in Fig.0). 



Thus, we define C s and Ch, s (resp. Cq and Ch,Q) as 
the angular correlation functions we obtain for the to- 
tal IGM signal and for the contribution due to patchy 
reionization through uncorrelated ionized bubbles within 
a uniform IGM when we only count in our model the bub- 
bles created by stellar radiation (resp. quasar radiation). 
In other words, we use the reionization history described 
in the previous sections (see Valageas & Silk 1999a) but 
to compute the CMB secondary anisotropies we only take 
into account the bubbles associated with either one of the 
two available sources of radiation (stars or quasars) . This 
allows us to compare the importance of stars and QSOs in 
our results (for the peculiar scenario of structure forma- 
tion we use). We show our results in Fig, 11. First, we note 
that we recover for both cases the main features described 



in Sect. 3. 3 for the total signal. Then, as expected, the 



comparison of Ch, s with Ch ; Q shows that the character- 
istic scale of the ionized bubbles associated with quasars 
is larger than for the HII regions produced by galaxies. 
However, the difference is not very large (note that the 
radius only scales as L 1 ' 3 , where L is the source luminos- 
ity, and one has to integrate over an extended population 
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of sources and over redshift). Moreover, we find that the 
total signals C s and Cq are very close and they do not 
exhibit different characteristic scales. Indeed, as we de- 



scribed in Sect. 3. 3 most of the power is provided by the 



small-scale matter density fluctuations of the IGM and by 
the large-scale cross correlations of ionized bubbles. Hence 
the typical size of the ionized bubbles cannot be seen in 
the shape of the angular correlation C{9). Besides, since 
quasars and galaxies are drawn from the same population 
of collapsed halos they have similar correlations hence the 
cross-correlations of their associated HII regions are rather 
close (see Valageas et al. 2000 for a detailed study of the 
correlation properties of these various objects). 




Fig. 12. The power-spectra S/ )S and S^q due to stars and 
quasars. The dashed curve Sj,h,s (resp. Si^q) correspond 
to the "homogeneous IGM" scenario for ionized bubbles 
created by stellar radiation (resp. quasars) within a uni- 
form medium. The dotted line shows the total IGM spec- 
trum Si (it is identical to the solid line Si in Fig.O). 



precisely, our model is based on the usual scenario where 
QSOs correspond to massive black holes located in the nu- 
clei of galaxies and powered by accretion (e.g., Rees 1984). 
Thus, an "exotic" model where quasars would not reside 
within massive virialized halos similar to galaxies might 
provide a different signature on the CMB. However, such 
a scenario is rather unlikely (e.g., in view of the energy 
requirements to power the quasars which favor large grav- 
itational potential wells) and the standard model has been 
shown to agree reasonably well with numerous observa- 
tions (e.g., the B-band luminosity functions and the X-ray 
emission, Valageas & Schaeffer 2000). As we have shown 
above, on small scales (I > 10 4 ) the signal comes from the 
fluctuations of the density field within the IGM and from 
galactic halos, while the pattern of reionization plays a 
minor role. Hence our results in this range do not strongly 
depend on the clustering properties of quasars. On the 
other hand, on larger scales the kinetic SZ effect probes 
the spatial correlations of QSOs and in this sense it be- 
comes more "model-dependent" . However, we can be rea- 
sonably confident in our results as our model has already 
been checked against observations of the QSO multiplicity 
functions (e.g., Valageas & Schaeffer 2000). Moreover, as 
shown in Fig. 8 in Valageas et al. (2000) we also recover the 
observed behaviour with redshift of the correlation length 
associated with quasars. This means that any model which 
satisfies the same observational constraints (up to z ^ 5) 
is likely to give analoguous results. Note that although the 
clustering properties of QSOs and galaxies as a whole are 
similar, since they are drawn from similar collapsed ha- 
los, the observed redshift-dependence of their correlation 
length is qualitatively different if one selects objects by a 
given luminosity threshold, due to the different behaviour 
of their mass-luminosity relations (see Valageas et al. 2000 
for a detailed discussion). 



We show in Fig.[l2| the power-spectra Si tS and Si t Q (as 
well as Si t h,a and Sui.Q for the "homogeneous" scenario) 
associated with stars and quasars, which also correspond 
to the correlations displayed in FigJO. We find again that 
the spectra Sfh >B and Sj,h,Q which directly probe the size 
of the HII regions exhibit two slightly different scales 
for quasars and stars, in agreement with FigJTP. Thus, 
Si,h,s peaks at I ~ 10 4,3 while Si^,q peaks at I ~ 10 38 . 
The wavenumber associated with quasar-driven bubbles is 
smaller than for stellar radiation since the size of the HII 
region is larger. However, we find again that this signa- 
ture is lost in the total power-spectra Si tS and Si t Q which 
are dominated at all scales by other processes (i.e. the 
correlations of the matter density field itself). Thus, ob- 
servations of the secondary anisotropies of the CMB are 
unlikely to provide strong constraints on the size of the 
ionized bubbles. Hence they cannot discriminate between 
both sources of radiation (stars versus quasars). 

Of course, this conclusion relies on the assumption 
that quasars are closely associated with galaxies. More 



3. 7. Dependence on cosmology 

Finally, in order to check whether our results strongly de- 
pend on the reionization history of the universe we also 
study the case of a standard CDM cosmology (SCDM): 
Q m = 1. We use the same model as in Valageas & Silk 
(1999a): fl h = 0.04, H„ = 60 km/s/Mpc and ct 8 = 0.5. 
The physical processes included in the model are the same 
as for the open cosmology and the galaxy and quasar mul- 
tiplicity functions are again normalized to observations at 
low z (Valageas & Schaeffer 1999; Valageas & Silk 1999a). 
The reionization redshift we get in this scenario is lower 
than previously: we now have z r j = 5.6. 

We show in Fig.O and in Fig.[l4| the angular correla- 
tion C{ff) and the power-spectrum Si we obtain for this 
critical density universe. First, we note that the ampli- 
tude of the secondary distortions of the CMB is smaller 
than for the previous cosmology. This is due to the smaller 
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Fig. 13. The angular two-point correlation function C(9) 
for a SCDM cosmology. The solid line labeled C shows 
the contribution from the IGM. The curve C ga i displays 
the contribution from galactic halos. The curves C u and 
Ch correspond to the "uncorrected bubbles" and "homo- 
geneous IGM" scenarios, as in Fig.||. 




log(l) 

Fig. 14. The power-spectrum Si of the secondary 
anisotropies for the SCDM cosmology. The solid curve la- 
beled Si (resp. S^gai) shows the contribution from the IGM 
(resp. from galactic halos). The curves Si tU and Si^ cor- 
respond to the "uncorrected bubbles" and "homogeneous 
IGM" scenarios, as in Fig|| The upper curve labeled "pri- 
mary" shows 1(1 + l)Ci/(2ir) of the primary anisotropies 
for the same SCDM model (see text). 



reionization redshift z r 
that: 



Indeed, the expression (Q) shows 



C{9) ~ Sz w(Zri) ~ Sz (1 + z ri ) 3/2 



(49) 



where Sz is the redshift interval during which reioniza- 
tion occurs. The smaller reionization redshift is due to the 
smaller variance as and to the faster growth with redshift 
of density fluctuations (which implies that for the same 
normalization at z — structure formation was less ad- 
vanced at high redshift for the SCDM cosmology). Next, 
we can check in the figures that we recover the same fea- 
tures for the secondary anisotropies as for the open uni- 
verse. In particular, we obtain the same shapes for all 



contributions to the correlation functions and the relative 
importance of each process remains unchanged. Hence all 
the points we discussed in the previous sections for the 
low-density cosmology remain valid. This is not surpris- 
ing since the basic astrophysical model we use is the same. 
However, it shows that our conclusions do not strongly 
depend on the details of the model (e.g., the reionization 
redshift). 

For these reasons, we expect that for low-density flat 
models, such as the 'cosmic concordance' ACDM model 
(Ostriker & Steinhardt 1995; Krauss & Turner 1995), the 
results should be similar to the OCDM model, except for 
a shift of the features to lower I (larger angular scales) due 
to the different angular geometry. 

3.8. Comparison with numerical simulations 

Here, we briefly compare our results with other avail- 
able studies. First, Benson et al. (2000) presented a semi- 
analytic model of galaxy formation to compute the kinetic 
SZ effect (within a ACDM cosmology). They roughly get 
the same amplitude Ci ~ 10 -13 at I ~ 10 4 and they also 
find that at low wavenumbers (/ Ss 10 3 ) most of the power 
is provided by density fluctuations and the clustering of 
ionization sources. However, because of numerical resolu- 
tion limitations they get a sharp drop at I ~ 2 x 10 4 while 
we obtained a plateau for Ci up to I ~ 10 6 . Note that 
our model is entirely analytic, although the scaling func- 



tion H(x) which enters the multiplicity functions in (17) 
is obtained from a fit to N-body simulations, so that we 
have no resolution limitations (we are simply limited by 
the approximations involved in our model). 

The secondary anisotropies produced by inhomoge- 
neous reionization which we study in this article have also 
been computed by means of numerical simulations in Br- 
uscoli et al. (2000), Gnedin & Jaffe (2000) and Springel 
et al. (2000), using different cosmologies and astrophysi- 
cal models. These authors find a broad maximum for the 
power-spectra Ci and Si of Ci ~ 10 -13 — 10 -12 around 
I ~ 10 4 . In particular, Springel et al. (2000) get a slowly 
decreasing plateau down to I ~ 500 (below this scale they 
are limited by finite box size effects while both other nu- 
merical studies are restricted to I *5 10 4 ). This behaviour 
agrees with our results (see the curve Si tU in Fig.g) since 
these authors use a simple toy model without including 
galaxy formation and radiative processes so that they miss 
the additional power due to the correlation of ionizing 
sources. Moreover, the drop we get for Ci at I ^ 300, where 
we recover a white noise power spectrum, is beyond the 
range of these numerical simulations. 

On small scales, Gnedin & Jaffe (2000) (with the high- 
est resolution) find a plateau at Ci ~ 10~ 13 which slowly 
decreases up to their resolution limit at / ~ 10 6 . This again 
roughly agrees with our predictions, although we rather 
obtain a slight increase of the total power with I in this 
range. Note that these numerical simulations are restricted 
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to z > 4 and these authors estimate the missing signal by 
a simple extrapolation (i.e. they multiply their output by 
a factor 1.25). However, as we discussed in Sect.[D] the 
redshift distribution of the kinetic SZ effect depends on 
the angular scale one considers. Thus, large wavenum- 
bers (I ~ 10 5 ), which probe high density fluctuations, 
are more sensitive to low z than large scales (I ~ 10 3 ), 
which probe the inhomogeneous pattern of reionization 
and where most of the signal comes from epochs close to 
the reionization redshift z^. This could explain the small 
difference between both predictions for the slope of this 
plateau. We can expect that with a higher resolution these 
authors would also recover a sharp cutoff at I <; 10 6 for 
Si (note that this drop is not readily apparent if one only 
computes the oscillatory spectrum Ci, see Fig.||). These 
two behaviours are recovered by Bruscoli et al. (2000) at 
l~2x 10 5 , but this smaller value for the location of the 
transition might be due to the finite numerical resolution. 
Since this scale is directly related to the size of virialized 
halos we can expect our result to be rather robust (in 
fact we would even expect some power at slightly smaller 
scales due to the collapse of baryons after they cool and to 
the substructures within halos, which would improve the 
agreement of our predictions with the results of Gnedin & 
Jaffe 2000). 

Bruscoli et al. (2000) also display the angular correla- 
tion function C(9). It reaches a plateau C(9) ~5x 10~ 12 
for 9 < 10~ 5 rad and it shows oscillations for 9 > 2 x 10 -4 
rad. Thus, the amplitude of the signal they get is larger 
than our predictions (as for Ci). This could be due in part 
to their higher reionization redshift, see (|49|). Moreover, 
Gnedin & Jaffe (2000) argue that those authors overesti- 
mate the SZ effect by a factor 3-10 because of the uncor- 
rected periodicity of the simulations. On the other hand, 
we obtain more large-scale power since in our model the 
cutoff of C{6) only appears for 9 <; 10~ 3 rad. On these 
large scales, secondary CMB anisotropies are generated 
by the cross-correlations of ionized bubbles. Hence this 
difference between both predictions may also be related to 
our smaller value for z T [ since in our case at reionization 
structure formation is more advanced and the correlation 
length of the matter density field (hence of the radiation 
sources) is larger. Moreover, these large scales are not ade- 
quately resolved by these numerical simulations, as shown 
by their results for C/ which are restricted to I <; 10 4 . 

Finally, we note that the independent study by Gnedin 
& Jaffe (2000) also finds that the signal is dominated by 
the contribution from high-density ionized regions rather 
than from the patchy pattern of reionization (for I > 10 4 ). 
This agrees with our results. Note that we find in addition 
that the inhomogeneous pattern of reionization plays an 
important role at larger scales (I ~ 10 3 ) but this is be- 
yond the range of these simulations. Thus, the agreement 
of our predictions with these various numerical studies, 
which use different cosmologies and astrophysical models, 



appears quite reasonable. Note that those numerical works 
do not include quasar formation models. 

4. Conclusion 

In this article, we have presented an analytic model (based 
on our previous work which described structure forma- 
tion processes and the reionization history of the uni- 
verse) which allows us to compute the secondary CMB 
anisotropies generated by the kinetic Sunyaev-Zel'dovich 
effect. This model includes a consistent description of 
galaxies, quasars and matter density fluctuations. 

We have found that the contribution due to patchy 
reionization is negligible except at very large scales {9 ^ 
10 -3 rad) and small wavenumbers (I Ss 10 3 ). Over this 
range, which corresponds to scales larger than the typi- 
cal size of HII regions, the signal actually comes from the 
cross-correlation of ionized bubbles, induced by the corre- 
lations of the rare radiation sources. On smaller scales, the 
IGM contribution is governed by the fluctuations of the 
matter density field itself. However, over this range the 
secondary anisotropies should be dominated by the con- 
tribution from galactic halos, which are characterized by 
smaller scales than the IGM (and larger densities). This 
leads to a cutoff of the power-spectrum l(l + l)Ci at a large 
wavenumber I ~ 10 6 . On the other hand, at low wavenum- 
bers I < 10 3 we recover a white noise power-spectrum. 
This very extended range of wavenumbers 10 3 < I < 10 6 
is close to the limitations of current numerical simulations. 
Thus, observations of these secondary CMB anisotropies 
should mainly probe the correlation properties of the un- 
derlying matter density field, through the correlations of 
the HII regions and the small-scale density fluctuations. 
We also found that the "local average" Si of the power- 
spectrum should be a more convenient tool than C;. 

Some comments are in order about the dctectabil- 
ity of the effects described in this paper. First, we no- 
tice that in the range 10 2 ^ I Js 10 4 the power pre- 
dicted by our model (relative to primary anisotropies) 
is comparable to the one found by Knox et al. (1998). 
Following their conclusions, we infer that this signal, if 
not taken into account correctly, might introduce a small 
bias in the determination of cosmological parameters from 
future experiments like MAP ( http://map.gsfc.nasa.gov|) 



and particularly Planck ( http://astro.estec.esa.nl/SA- 
[general/Projects/Planck/ ) . Second, although the range of 



I where our model produces most of the power (10 4 Ss 
I & 10 6 ) is likely to be out of reach for MAP and 
Planck, future mm-wavelcngth interferometers, such as 



ALMA (http://www.mma.nrao.edu) may have the right 
sensitivity (~ 2/iK rms for a 1' beam in one hour) and 
the right resolution (< 2') to be able to measure such a 
signal. Indeed, although the amplitude of the secondary 
anisotropies which we obtain is somewhat lower than the 
sensitivity of ALMA, a larger normalization of the power- 
spectrum (for the SCDM case we used a$ — 0.5 while 
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a COBE normalization would give <j% ~ 1.4) or a larger 
reionization redshift would push the signal into the range 
of detectability. 

We noticed that the redshift distribution of the con- 
tributions to these secondary CMB anisotropies is rather 
broad. Thus, for the angular correlation from the IGM we 
get 7.5 < z < 10 with a sharp cutoff at the reioniza- 
tion redshift z r ; = 6.8, when the "patchy pattern" of hy- 
drogen ionization disappears. However, some small-scale 
anisotropies are still produced at lower redshifts. The red- 
shift distributions of the contributions from galactic halos 
are even broader, we typically get < z < 7, and show no 
strong feature at z r {. Since the total signal should be dom- 
inated by the contribution from these collapsed objects for 
a large range of wavenumbers (I ^ 10 4 ) this implies that 
one should not assume that most of the secondary CMB 
anisotropies arc generated during a small redshift interval 
Sz around the reionization redshift z T [. 

Next, as expected we have found that within our sce- 
nario ionized bubbles produced by quasars are larger than 
those built by galaxies. This implies that the "patchy 
patterns" of the HII regions associated with QSOs and 
stars are different. However, since the total signal is dom- 
inated by the correlations of the matter density field, and 
not by the size of the ionized bubbles, it is similar for 
both radiation sources (which also have similar correlation 
properties). Hence, unfortunately one cannot distinguish 
a quasar-driven reionization process from a galaxy-driven 
reionization history, using the CMB anisotropies. 

Finally, we have checked that our predictions apply 
both for an open cosmology and for a critical density uni- 
verse. Thus, our conclusions do not depend on the cos- 
mological scenario and can be extended to low-density 
flat models. However, the amplitude of the anisotropies 
is larger for the low-density universe because of the higher 
reionization redshift. 
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